PRA-HEP/97-7 



Jet rates in deep inelastic scattering: wherefrom comes the sensitivity to a^? 

Jin Chyla and Jin Rames 
Institute of Physics, Academy of Sciences of the Czech Repubhc, Prague [| 
. . . , Abstract 

. 

a^ ■ 

0^ ' For theoretically consistent determination of ag from jet rates in deep inelastic scattering the 

• dependence on Us of parton distribution functions is in principle as important as that of hard 

scattering cross-sections. For the kinematical region accessible at HERA we investigate in detail 
^ I numerical importance of these two sources of the dependence of jet rates. 

^ ; 1 Introduction 
T— I ■ 

^ I One of the problems of quantitative determination of the running of is related to the fact that it 

' usually requires combining results of different experiments in different kinematical regions. Recently 

. the HI Q and ZEUS ||2| Collaborations have reported evidence for the running of as(^) obtained 

i from the measurement of jet rates in deep inelastic scattering (DIS) via the quantity 

O . 

0> ■ 7? ^^2^ _ Q-2+l(Q^) . ^ 

^ . 0^1+1 (Q^) + 0-2+1 (Q^) 

^ . where cjfc+i denotes the cross-section for the production of k hard and one proton remnant jets. 

D I There are now several NLO Monte-Carlo generators ^, ||, suitable for analyses of jet production 

' in DIS. In this paper we shall concentrate on the HI analysis Q which, using PROJET 4.1 |Q with 

^> ■ the JADE jet algorithm and yc = 0.02, obtained the following result for as(M^, MS) 

V . 

U ■ as{Mz, MS) = 0.123 ± 0.012(stat.) ± 0.008(syst.). (2) 

In JADE cluster algorithm, the jet resolution parameter yc is defined as yc = {pi -\-Pj)'^ /W^ , where 
Pi,Pj are parton momenta and W'^ = Q'^{\ — x)/x stands for the square of the 7p CMS energy. 
The result @j has been extracted from the Q^-dependence of i?2+i(Q^), displayed in Fig. 2 of Q, 
using only the two highest data points. Within the error bars the value (||) is consistent with 
the world average, but its accuracy is insufficient to draw any firmer conclusions. In the procedure 
adopted in as(/x, MS) (or, equivalently, Aj^) was considered as a free parameter in the hard 
scattering cross-sections, but not in the parton distribution functions (PDF), for which the MRSH 
set was used. However, as each set of PDF has a particular value of Aj^ built in, one must check 
for the consistency between this input Aj^ and the output one, obtained from comparison of (|^) 
with data. Clearly, in a consistent determination of a^, Aj^ must be varied simultaneously in PDF 
and parton level hard scattering cross-sections. In practical applications to physical quantities 
involving beside the hard scattering cross-sections also parton distribution and/or fragmentation 
functions it is often useful to know wherefrom comes most of the sensitivity to Ug ■ The purpose of 
our paper is to investigate this question in detail for the jet rates (|l|). 
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Throughout this paper we stay in the conventional MS renormahzation scheme (RS) of the 
couplant a = as{n)/7r and omit therefore the specification "MS" in as weh as in as(/i, MS). 
The jet cross-sections ak+i in (|l|) are given as convolutions 



, Jo 



x,ycA) 



(3) 



of the parton level cross-sections C^+i^j and PDF fi{x, M), evaluated at the factorization scale M. 
The sum runs over all parton species i. In perturbative QCD Ck+i,i are given as expansions in the 
couplant a{p,/A), taken at the hard scattering scattering scale |^ /i. To the NLO we have 



C2+i,i{Q, M, X, Vc, A) = a(^/A) 4^i,i(Q^ x, y^) + a{f^ / A^^^^/Q , M, x, y^) 



(2) 



Ci+i,i(Q,M,x,ye,A) = cfi,,(g,x) + a(^/A)clYi,i(Q,M,x,y,). (4) 

In standard global analyses of hard scattering processes §], A is fitted together with a set 
of parameters af \ describing distribution functions p^^\x,M) of parton species j at some initial 
scale Mq, usually in the form 



p^^\x, Mo) = (1 - x)4'' (l + a'i^V^ + a!-px + afx^'^) 



(5) 



Writing the derivative d(T2+i/dlnA as 
da2+i(Q2,A) 



E 



dlnA 

1 

dx 



dlnA 



{M)Y: f'dx \- [' ^My,M,A)P^p{z)a{M)c^^l^^{x)+ 
~7 Jo . Jo y 



bfi{x,M,A)a{M)4l,.{x 



(6) 



where z = x/y and P-^j (z) are the LO branching functions 0, we see that the leading order term 
of dcj2+i/dlnA gets contributions from the variation of A in both the PDF fi and hard scattering 
cross-sections C2+i,i- The two terms in the brackets of @ are of the same order and their relative 
importance thus is basically a LO effect. 

For the quantity i?2+i((5^) the situation is less obvious. PDF appear in both the numerator and 
denominator of (|l|) and so some cancellations might occur, while the hard scattering cross-sections 
start as 0{a) for C2+i,i and as 0(1) for Ci+i^j. Nevertheless a detailed analysis |1C] shows that 
also for the ratio i?2+i varying A in PDF generates terms which are of the same order as those 
resulting from the variation of A in the hard scattering cross-sections C^+i^i- 



2 Parton distribution functions for arbitrary A 

To assess the potential of (||) for a precise determination of Os we wish to separate the question of 
its sensitivity to as from the sensitivity to parameters describing the initial condition on the PDF. 

^Although in general fj, ^ M, we shall follow the usual practice and set fi — M. 

^In (^) we have for brevity suppressed the dependence of c^^^ ^ on yc and written a{M/A) simply as a{M). 
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Moreover, in order to investigate the numerical importance of varying A in PDF we have to know 
what to do with the parameters Mo,ap^ in (^) when A is varied. We cannot keep them fixed as 
this would contradict the idea of factorization. To see why, let us first consider the couplant as well 
as the evolution equation for the nonsinglet quark distribution function q^s{x,M), at the LO. In 
terms of conventional moments we have 

^^^^^^^^ = aiM/A),.sin,M)pM{n), = -<s^(n)/6, (7) 

whence 

QNs(n, M, A) = [a(M/A)]"" , a(M/A) = ^^^^^^y (8) 

In the above expressions pj^g (n, M) are moments of LO nonsinglet branching function and An are 
unique dimensionless constants, describing the nonperturbative properties of the nucleon, intro- 



duced in 1 11]. Via (|^) they determine the asymptotic behavior of q-!<is{n, M, A) for M — > oo and 
thus provide alternative way of specifying the initial conditions on the solutions of The sepa- 
ration in ^ of q^^si'iT': M, A) into two parts - one calculable in perturbation theory and the other 
incorporating all the nonperturbative effects - is the very essence of the factorization idea. Note 
that in (^) the dependence on A is simply a reflection of its dependence on M. Knowing the latter, 
we know the former. Eq. (^) implies that we can write down the equations for the A-dependence 
of PDF which are very similar to standard evolution equations, describing their dependence on the 
factorization scale M. The situation is basically the same as for the running couplant a(M/A) for 
which do(M/A)/dlnA = — da(M/A)/dlnM. In the case of PDF the only difference stems from 
the necessity to properly specify the initial conditions on the solution of For finite initial Mq 
implies 



qt^s{n,M,A) = qj^s{n,Mo,A) 

where 



a{M/A) 
a(Mo/A) 



gNs("',Mo,A)exp(-(i„s), (9) 



is the so called "evolution distance" |Q. The second equality in (|T0|) holds exactly at the leading 
order and approximately at higher ones. Rewriting (^) again in the form (P) 



QNS (n,M, A) 



q^s{n,MQ, A) 
{a{Mo/A)) 



d„ 



{a{M/A)f" (11) 



we see that the ratio of q^sin-, Mq, A) and (a(Mo/A))'^" equals An and must therefore be both Mq 
and A independent. This, however, is possible only if the initial moments qjsisin', Mq, A) do, as 
indicated, depend beside Mq on A as well! Note that considering q^sin-^ M, A) as a function of 
s, formula ( p!I| ) implies that entire dependence of q]<fs{n, M, A) on A comes solely from the term 
ln(M/A) = l/a(M/A) in (0)! The same happens in the realistic case of coupled quark and gluon 
evolution equations ||T^. If we wish to investigate the dependence of PDF solely on A , keeping 
fixed all other parameters, specifying the initial conditions on PDF, we cannot fix initial conditions 
at some finite Mq, (i.e. the parameters p^^\x, Mq) in (||)). It is the constants An that must be kept 
fixed instead! On the other hand, in global analyses, like 0, ^, all parameters, including A , are 
varied simultaneously and fitted to data. In these circumstances it is then not straightforward to 
determine the sensitivity of a given physical quantity to itself. 
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The preceding paragraph also suggests a simple procedure, which makes use of some of the 
available parameterizations and generalizes them to arbitrary A , keeping the constants An fixed. 

From all the available parameterizations of PDF those given by analytic expressions of the coef- 

(i) 

ficients af on a general factorization scale M via the variable s are particularly suitable for this 
purpose. In our studies we took several of the GRV parameterizations, both LO and NLO ones, 
defined in pi. At the LO the recipe for construction of PDF for arbitrary A is simple: 



• Use any of these parameterizations with fitted value of A fn . 

• Keep A in ln(Mo/A) fixed at the value Aju. 

• Vary A in ln(M/A). 

• Use the original parameterization with the redefined s = s{A). 

This construction is exact at the LO. At the NLO the couplant a(M/A) is no longer given simply as 
l/ln(M/A) and consequently the second equality in (|lO|) holds only approximately. This, together 
with the additional term appearing in the NLO expression for q-f^sin, M, A) 



qmin,M,A,An) = An 



a{M/A) 1 , 



1 + ca{M/A) 



(l + ca(M/A))^NsW/^'=, (12) 



destroy simple dependence of ^ns on M, Mq and A via the variable s. Using the above recipe for 
the NLO GRV parameterizations provides therefore merely an approximate solution of our task, 
but one that is sufficient for the purposes of determining the relative importance of varying A in 
PDF and hard scattering cross-sections as this relative importance is basically a LO effect. 

As an illustration of our procedure we plot in Fig. 1, for x > 10~^ and two values of Q^, 
the dependence of valence and sea u-quark and gluon distribution functions on A. The curves 
correspond to the LO GRV parameterization of p!3[| . We see that this dependence decreases with 
increasing Q^, is most pronounced for the gluon distribution function and almost irrelevant for the 
valence quark one. These features are qualitatively the same as those obtained in refs. |14, 15 1, 



which contain results of global fits performed for several fixed values of A. Quantitatively, however, 
the sensitivity to the variation of A , observed in these papers, is markedly weaker than that 
displayed in Fig. 1. This is due to the fact that our A-dependent PDF are not constructed to 
describe the hard scattering data for all values of A (they do so by definition only for A = Afn) 
but in order to facilitate studies of the dependence of physical quantities solely on A , with the 
appropriate parameters specifying the initial conditions kept fixed. On the other hand, in global 
analyses of |l^ fitting the data for different values of A leads to different values of boundary 
condition parameters of '' at Mq, which partially compensates the variation of A. 



3 Numerical results 

All the results reported below were obtained using the PRO JET 4.1 generator |^. It is true that 

the JADE jet algorithm used in this program has certain theoretical shortcomings, but as the HI 

Collaboration used it, we did the same. Although MC generator MEPJET Q is to be preferred 

on theoretical grounds, the main features of our results are unlikely to change because the A- 

dependence of R2+1 is basically a LO effect. To quantify the dependence of R2+1 on A separately 

in hard scattering cross-sections and PDF we consider it as a function of , yc and two independent 

A-parameters, denoted as A'^^^ and Ap*^^ respectively. The results are then studied as a function of 
^hsc ^ ^pdf ^ q2 
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Figure 1: The dependence of xuy{x , Q"^) ^ xusea.{x , Q"^) and xG{x,Q'^) of the value AP*^^ according 
to the recipe described in the text and for two values of Q^. At low x, the curves correspond from 
above to the values of AP'^^ , subsequently, 0.1,0.2,0.3,0.4 and 0.5. 
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Figure 2: The ratio ii2+i(<5, yc, AP*^^ ,A'^'^'^ ) as a function of AP*^^ for fixed A^^*^ and a series 
of values (solid lines), or as a function of A'^'^'^ for fixed AP'^^ and the same set of values 
(dotted lines). The dashed lines correspond to simultaneous variation of AP'^^ = A'^^'^ . Solid curves 
correspond, from below, to fixed A^^'^ = 0.1, 0.2, 0.3, 0.4, 0.5 and the dotted to the same fixed values 
of AP'^f . In all plots yc equals 0.01. 

• 12 values of (equidistant in InQ^) between 5 and 10^ GeV^, 

• 5 values of AP'^^ or A'""'^ , equal to 0.1,0.2,0.3,0.4 and 0.5 GeV, 

• 3 values of = 0.01, 0.02, 0.04. 

3.1 No cuts on jets 

First we analyse the case with no cuts on the final state jets and then discuss the changes caused by 
the imposition of cuts as specified in |]l|. Fig. 2 shows a series of plots corresponding to yc = 0.01 
and displaying the dependence of i?2+? A^'^'^ for fixed AP'^^ (dotted curves) or on AP'^^ for fixed 
A'^**'^ (solid curves). Each plot corresponds to one of the 6 selected values of Q^. We see that for 
Q2 below about 40 GeV^ the solid curves are steeper that the dotted ones, while above 40 GeV 
the situation is reversed. This means that below 40 GeV^ varying A in PDF is more important 
than varying it in the hard scattering cross-sections for Q^, while the opposite holds above 40 
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Figure 3: (a-c) Q'^-dependence of the ratio i?2+i(Q, yc, AP'^f , A^"^ ) for fixed AP^^^ = 0.1,0.3,0.5 
GeV and five values of A^'^'^ = 0.1,0.2,0.3,0.4,0,5 GeV. (d-f) The role of Ap^^ and A^^^^ is 
reversed. No cuts were applied in evaluating R2+1 and yc = 0.01. In all plots the curves are 
ordered from below according to increasing A^^'^ (in a-c) or AP'^^ (in d-f). 

GeV^. The same information as in Fig. 2 is displayed in a different manner in Fig. 3, where the 
dependence of R2+1 is plotted for various combinations of AP'^^ and A'^^'^ . The shape of the 

curves in Fig. 3 is a result of two opposite effects: as increases the phase space available for 

produced jets increases as well, but at the same time the value of the running as{Q/A) decreases. 

At low Q^, on the other hand, phase space shrinks, but as{Q/A) grows. 

In Fig. 3 the relative importance of varying Ap*^^ and A'^'^'^ is reflected in a bigger spread of the 

five curves in Figs. 3d)-f) compared to those in Fig. 3a)-c) for < 40 GeV^ and smaller spread 

above 40 GeV^. Note also that for fixed AP'^^ the sensitivity to the variation of A'^'^'^ is about the 

same at all Q^. The shape of curves in Fig. 2 suggests that the ratio 

,.,^2 , _ fi2+l(Q^yc,APdf =g,Ahsc =r) 

= i^2.i(Q^y.,Ap^^ =r,A^- =sy ^''^ 

is approximately a linear function of s for any fixed Q'^,yc,r. Considered as a function of s for fixed 
r, (|l^ quantifies the relative importance of varying A in PDF and hard scattering cross-sections. 
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Figure 4: (a-c) The quantity W{Q'^,yc,r) as a function of for five values of r (given in GeV) 
and Uc = 0.01, 0.02, 0.04. All curves correspond to the case of no cuts. 

To summarize the results of Fig. 2 we fitted V{Q'^,yc, r, s) by a linear function of s and defined 

WiQ ,yc,r) = — , (14) 

which, by construction, is s-independent function of Q'^,yc and r. Positive W means that the 
variation of A in the PDF is more important than that in the hard scattering cross-sections, while 
for negative W the situation is opposite. The Q ^-dependence of W{Q'^,yc, r) is plotted in Fig. 4a-c 
for three values of yc = 0.01, 0.02, 0.04 and five values of r = 0.1, 0.2, 0.3, 0.4, 0.5 GeV. In Fig. 5a we 
compare the Q2 dependence of W{Q'^,yc,r = 0.2 GeV) for different values of yc = 0.01,0.02,0.04. 
Similar plots can be drawn for other values of r as well. From Figs. 2-5a we conclude that at the 
NLO: 

• For Q2 below 40 GeV^ the sensitivity of (|I|) to A comes dominantly from PDF, while above 
40 GeV^ the situation rapidly changes and most of this sensitivity comes from hard scattering 
cross-section. 

• Above ~ 10^ GeV^ the sensitivity to A in PDF becomes negligible. 

• The preceding conclusions depend only weakly on y^.. 

At the LO the general features are the same as those displayed in Figs. 2-5a and we therefore 
merely summarize them in Fig. 5b. Comparing Figs. 5a and 5b we conclude that at the LO the 
relative importance of varying AP'^^ with respect to A'^^'^ is bigger than at the NLO. 
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Figure 5: The Q^-dependence of W{Q'^,yc,r = 0.2) for three different values of yc at the NLO (a) 
and LO (b). All curves correspond to the case of no cuts. 
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Figure 6: The same as in Fig. 3 but for the HI acceptance cuts. 



3.2 HI cuts 

In IQ] only events satisfying certain acceptance cuts, most notably the cut on polar angle of outgoing 
jets in HERA laboratory frame, 10° < Olf^ < 145°, were selected for the analysis. In Fig. 6 
plots analogous to those of Fig. 3, but incorporating these cuts, are displayed. We see that the 
experimental cuts affect mostly the low region, where beside lowering the values of i22+i(Q^) 
they also strongly suppress the sensitivity of i?2+i(Q^) to the variation of A in the PDF. This effect 
can be understood as follows. By imposing the mentioned angular cuts, one removes events with 
jets that fly into a large part of the backward (with respect to the proton direction) lab frame 
hemisphere. We recall that in direct photon interactions (in contrast to the resolved ones), i.e. 
also DIS, a large fraction of jets populates the backward hemisphere in the 7p CMS. Because the 
transverse boost between the 7p CMS and HERA lab frame is proportional to Q^, the cut on the 
jet polar angle in the lab system is more effective for small than for high Q^. Hence the drastic 
suppression of -R2+i(Q^) at low and only moderate at high Q^. Furthermore, as low means in 
average also small x, and most of the sensitivity of i?2+i(Q'^) to A in PDF comes from the small x 
region, the HI cuts will significantly reduce it mainly for low Q^. 

On the basis of Fig. 6 we conclude that in the region > 100 GeV^, used in Q for the 
extraction of A from the measured ii2+i(Q^), the variation of A in the PDF could be neglected 
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with respect to the variation of A in the hard scattering cross-sections. 



4 Conclusions 

In this paper we have constructed parameterizations of PDF corresponding to solutions of the 
evolution equations with fixed suitably defined boundary conditions and variable A . These pa- 
rameterization were then used to investigate, for the quantity -R2+I) the numerical importance of 
varying A in PDF. We concluded that: 

• If no cuts are applied on produced jets the sensitivity of (||) to variation of A in PDF is bigger 
at the LO than at the NLO. 

• At moderate Q^, roughly below 40 GeV^, the sensitivity of (|^) to comes dominantly from 
PDF, while at higher the hard scattering cross-sections take rapidly over as the main 
source of the A-dependence of (||). 

• The preceding conclusions are only weakly dependent on Uc- 

• Including the cuts applied by the HI Collaboration leads to strong suppression of the sensi- 
tivity to Us in the moderate region. 

It is clear that a less restrictive treatment of the jets in the region between, say, 10 and 100 
GeV^ could provide much better possibility for a precise determination of a^. This region is also 
very interesting from another point of view, and namely the question of the transition between the 
dynamics of deep inelastic scattering and photoproduction. The data in this region could bring 



new information on the structure of the virtual photon [16| and/or lead to new phenomena. 
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